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I. SEVERAL VORTEX EI LAMENTS ■ 



In the airfoil theory of Prnndtl (reference l) the 
trin^ is replaced by a lifting vortex filament whoso circu- 
lation v.iries over the span. By this aethod -he " flrat 
problem of airfoil theor-'," namely, for a ^lven lift dis- 
tribution to dot ermine the shnpc of the airfoil, was 
solved. The inverse "socond pro'blom," nnmoly, for n -Sivon 
rin^ to determine the lift di at ri "but ion , vas then solved 
by Bets (reference 2), the computation being: simpler for 
small aspect ratios than for lar^e ones. For the latter, 
an approximate solution was obtained "by Trefftz (refer- 
ence 3). The answer vras thus found to the most important 
practical question, namely, the manner in which the rrin^ 
forces are distributed alon^ the span. 

The chord^ise distribution theory was simply taken 
over fron the theory of the infinite win^. The Ackermann 
formulas, published by Birn"baum (reference 4), in which 
the infinite \rin«; was replaced "by a plane vortex sheet — 
on account of their linearized form permit also applica- 
tion to the finite win^ and this ax>pli cation vied' carried 
out by Blenk (reference 5) for the rectangular win^. 
Since in this work a series expansion in "b/twas. used, 
the computation converges only for lar^e aspect ratios. 
In the present paper a useful approximate solution will be 
found also for win<$s with large chord - i.e., small aspect 
ratio . 

Another method of investigating the lift distriou-. 
tlon alons the two dimensions (span and chord) was found 
"by Prandtl, (reference 6)' in his use of the acceleration 
potential." This method assumes, horever , -'that the poten- 
tial is known for a suitable number of source distrihu- 



* n Uher die Auf triebsverteiluns des einfacaen Rechteck- 
flu^els uber die Tlefe." Zaitschrif t fur an^ewandte 
Mathematik und Ilechanik, vol. 19, no. 5, Oct. 1939, 
pp. 257-270. 



2 



HACA Technical Iieatfrandun JTo.- 983 



tions over the horizontal projection.- of the win^» The 

first application was made by Kinner (reference 7) in hie 

work on the win^ with circular plan form, since these 
functions are obtainable for the circle, The nethod ap- 
pears, however, for the present to offer no promise for 
the rectangular wing, since no expansion of the potenttal 

into a series of known functions la known for the rec- 
tangle. Bar this reason the computation In the present 
paper will still be conducted b7 the vortex-shset method. 

For accurate investigation of the lift di stri"bution , 

the wing nust be replaced by a vortex sheet. A qood idea 
of the distribution can still "be obtained if the win*; is 
represented "b^ a finite nunber of discrete vortex filaments, 
and the necessary amount of computation is there-by reduced 
considerably as compared with the continuous circulation 
distribution. This is because in the case of t^e vortex 
sheet, the condition that the component of the induced ve- 
locity at right angles to the win? should be equal to that 
due" to the flow, ?iven rise to en integral equation. Bor 
individual vortex filaments, however, this flow condition 
need be satisfied exactly only at single points, so that 
only a system of linear equations is obtained. Figure 1 
shows such a vortex system, for which the computation was 
carried out. In order that the results obtained from 
usinfl; only a few vortices, or euere a single one. be as 
accurate as possible, the distance of the first vortex 
filament from the leading ed?e is taken to be a/4. It is 
knom fron previous wori that the circulation in the neigh- 
borhood of the leading ed?e increases as l/-/xT the fore- 
most, strongest vortex aiich ?ives the circulation contribu- 
tion fron tho leading ed?e to the foremost points consid- 
ered, then lies exactly at the center of pressure of the 
forward lift portion because the center of gravity of y = 
o/i/x lieu at s = x/3. The points at which the total ve- 
locity at ri?ht angles to tho win? is nado to vanish, lie 
in the oonter betwoen two vortex lines and at x = t a/4. 
Tho win<5 is a plane rectangular pinto of zoro thickness 
with chord t = na for n vortices. Tho notation is in- 
dfcated in figure 1. The coordinates of_tho po*nt A aro 
x* and y*. The velocity at ri^ht angles to J no xy 
plane, inducel by the bound and trailing vortices at the 
point A, is thon ?iven by the Biot-Savart law: 
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The fir-t integral, which arises from the bound vorticos, 
■^Ives after integration "by parts: 



|" (y - y*) r^Cy) i + 



+ £>/? 

V2 



+T)/2 



-/ 

-o/2 



(r - y«) 



dr ± ( y ) 



^ y Xi + (y - y*) 



dy 
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The first expression vanishes since, the circulation at the 

tips Lust be zero; r (± ^ ) = 0. There is thus obtained 
"A: * 2/ 

+t>/2 



1 J dy y-y* ^ Sj _ ' 



dy (la) 



-h/2 
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Since (y) decreases from the center of the plate to- 

ap. (y ) 

ward the tips, — |^ — < 0, and hence vt± > 0 . Iron tho 

condition WA = V sin a, there are obtained for n points 
m oquations xr± = = . . . = a 7 sin a. It is thus 

possible either toassune^oeano spanwlso circulation 

distribution - for exanplo, the elliptic for all n = n 
vortices — or set up a series expansion rrith r undeter- 
mined coefficients for n * m/r vortices. 

An example by the second method will ffrst "oo comput- 
ed. For this purpose, the following transformation of 

coordinates is made: 

y = -|cos q>i y* = |cos cp*j Xi = x* - x A = | fij 

so that - | 5 y 5 ^ corresponds to tt > cp > 0. Equation 
(la) then becomes: 



2ttd \ J COB cp - cos cp* 



( J*'* + 'cos cp - cob cp*T*\ dF, (cp) 

V- + — * S ) - dep (2) 

' dco 

For each T* (cp), a triqononet ric series that contains on- 
ly thesin (2v + 1 ) cp terns ras assur.ed, since the rela- 
tions are assumed symmetric with respect to the win 1 ; cen- 
ter: Tj_(cp) = cin cp (1 + a^^sin cp + a^ 8 ^ sin 3 cp ) . 

Thus for each vortex f ilnnent , there are three undetermined 

tx) (a) 
coefficients Ti , aj, * , and aj. . Be the- h\ *e 

dTi(cp) n , li) \ 

— =~ — = Ti (cos cp + 2aj; ' sincp coscp 

+ a^ 3 ) cos cp sin.Sp + 3a^ 8 ^ sincp cos 3cp) 
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Substituting this expression in equation (2), there Is ob- 
tained: 

W A = - Frk f 1 r i l^i J - 71 + a i ( ° f i (6 i« cos ^ - 

+ ai (s5 f a - (8 it cos qp*)} (3) 

where the functions t x and f a are made up .of integrals 
which may be evaluated "by elementary methods. ' The inte- 
gral J, which i"s also a function of8^ and cos cp* is, 

in any particular case to be determined- ty graphical or 
numerical methods. The integral is 

TT 



j 3 f cos qp dqp^ H 

J cosqp -coscp * 



* \3 



± + (COS Cp - COB Cp ) 



TT 



(cos cp - cos qp) , , . 

■ — — — coa cp dcp + tt 6. 

« 8 . f & a , STI 1 

°± + «/ 6 i + (cos cp - cos cp ; 



for 8jl and cos* > 0. Through this trr.nef ornation the 
singularity at cp -» cp* h.is bean removed. Since for each 

vortex line there are three undetermined coefficients, the 
flow condition can be satisfied for each set of three 
points "between two lines and for three joints at x = 
t — a/4. On account of tie symmetry 3n^ different points 
□ ay be chosen on a half win*; and for the corresponding 
points, symmetrical with respect to the center line of the 
plate, the condition WA »r Tf sin a is then automatically 

satisfied. Altogether, therefore, the flow condition is 
acouratelp satisfied for 6n points or, in case one of 
each set of three' points lies on the centor line, for 5n 
points. The entire computation is based on the expecta- 
tion that the condition w^ * 7 sin a will be, on the 
average, satisfied at least approxinat ely , also at other 
points of the surface, and that the singular behavior of 
Vjl along each of the lifting vortices will not havo too 
^j^bjijA ah effect on the approximate computation of the cir- 
oul'atl although norodynaiiically this can only be jus- 
tified':; oonB&aering the r>late as ro^lqced by several 
winga lying one behind the-other, each represented by a 
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vortex filanent. The choice of the number of vortex fil- 
aments n is, for practical reasons , rostrictcd sinco- 
while the number of points cocsidared in.oreasos only lin- 
early& with n 3 the required conputation Fork of solving 
tho system of 3n equations increanoe at a creator rate. 

The numerical commutation was carried out for the fol- 
lowing case: n = 4, q> = 30° (150°), 50° (120°), and 90° 

(center line) with b = 4a, corresponding to an aspect 
ratio of A = b/t = 1. (See fiq. 2.) The integral J 
(J (+ 6^) = J (- 8 4 )). was determined for the four values 
1/4, 3/4, 5/4, 7/4 which are assumed "by 84 and, on ac- 
count of the symmetry, for only three values of COB cp*. 
For qp* = 90°, an elliptic integral of the second kind is 
obtained for J. Pror. the functions f 1 (^1, cos qp ) and 
fj (Sj, cos "qp 5 }, the coefficients were obtained for a Bea- 
ten of 12 equations which was solvod by the usual elimina- 
tion process with the conputation machinesinco tho systom 
could rot ho solvod "by iteration. As tho conputation was 
carried out to only five decimal placos, it was afterwards 
found to be of .insufficient accuracy for tho determination 
of tho last throe unknowns; tho circulation of tho rearmost 
vortex filanont, therefore , coals only 00 estimate? "by ex- 
trapolation. For the romainin^ circulations, thore w n .s 
obtained: 

I\ = +0.737 0 bV sin a T a = +0.116 3 hV sina T 3 = *K).058 43 bv sina 
a^ 1 ) = -0.136 3 s= -*0.564 a &^ = +0.245 8 

-K).005 ? a g a ) ^ +0.001 a a£ a ) = -0.063.J 

Those circulation distributions are shown on figure 3. 
Integrating ovor the span: 



TT 



P 7 I y\to) sin (p d(P = p 7 I r t (f + § a^ 1 ) - A a;<">) 

0 

these mav ho considered "by tho Kutta-Joukowsky thooron as 

tho lift contributions of the individual win^ strips (alon^ 
the chord), Tho lift is then distributod as follows: 



NACA Technical Memorandum Ho. 963 



7 



Iron the leading ed^o to 3t/l6 : Aj -0.512, pbtV 8 sin a 
Prom 3t/l6 -to 7t/l6 : A IX = 0,133 B ■ 
From 7t/l6 to llt/l6 : A III = 0 055s " 
From llt/16 to trailing ed.j;e A I7 ~ 0.02 

The total lift is A = EpA^ = 0.72s p "b t 7 s sin a and 

' to . 

the lift coefficient c a = A / - "b t 7 = 1.44, sin a. The 

moment about the leading ed^e is obtained as the sun of 
the moments of the several strips: M = t/16 (AI + 5Ajj + 

9Aj j j + 'I3A17) ond the moment coefficient cm = 0.24s sin a. 

The -.position of the center, of .pressure is' obtained from 

9 = — 1 t = 0.16, t, 'vhoro s is tho distance of the cen- 
c a 

ter of pressure from the leading edge. Jinally, for the 

TT 

/ T ± (cp) sin q> dep 

forward three vortices the factor 7*4 3 - n r.-* m — 

can lac determined: = 0.80 9 , v a = 0.74,. v 3 = 0.73,. 

Sinco for all vortices this factor is approximately 
'equal to tt/4 = 0.78 64t it appoars justifiable — at least, 

for deep wines, that is, small aspect ratios A = b/2 — to 
assvime initially an elliptic sp&nwiso lift distribution and 
so considerably simplify the computation. The fact that 
this assumption . according to the above cor.puted example, 
is not quite applicable to the roar vortices, is of no great 
importance on account of the strong rearward drop in the 
circulation, 

Since for elliptic di st ribution there is only one un- 
determined coefficient for oach sorter lino namely , the 

circulation T*_ in the center of the span, the flow .condi- 
tion = V sin a can also "bo satisfied at only one point 

(and the point symmetrical with resgnct to the center line) 
■between each two lines and at - = 15t/l6 t For this rea- 
son, the assumad points are ta^-or. on tie center line and in 
the center between each, two sucr.codin !T . lines, and the last 
at x = 15t/l6. Te then have y* = 0, !*■= i (zjc+Xfc+i ) 1 
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, — t at 

for the rearmo st point x*= 15t/16 and x^ = *2~n' * g^* 

• • + ^Zn ~ ^ * where n is the number of vortices. Substi— 

8 r 4 y 



, „ / /y n b dri (y) 

tatm< ri (y) = m/i - (^-;. -^-= 



\("b/2 )/ 



into eguation (la) , there is obtained for the induced veloc- 
ity at the point considered: 



"A ■ 4» zl Fi A -7—- a I : 



1 dy 



+ 0/2 



= ~~rs h x i In - + — / / - 1 dy ^ 

ttd l 2 xj. J / / y »■ 



The elliptic socond integral is roducod to the normal form 
or the substitution y = b/2 costy, sothat there is ob- 
tained: 

where E is tha complete elliptic integral of the second 

tt/2 



(k, |j = y /l -k a sin a v|/ dV with 

o 

This function is tabulated, for example, 



kind E (k, -J = y ^ 1 ~k sin v|/ d\|/ with modulus k « 

o 

a/g 

fxf + (h/2) J 

in Jahnko-Em&e : Funktionentcf eln . it has the following 

limiuin^ value for k = sin <p, 0 < k_£ 1, 0 < cp < tt/2, 

and in the °.oovo equation for » > | xjj > 0: ~> E ^cp, HJ^i 
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The expression in parontheses in equation (4) gives 
the coefficients for the linear nonhomogeneous system of 
equations for tbe n' unknowns r x , r a i • • Tn m Cutting 2b 

on the right side, it reads 2b V sin a for all equations, 
on account of the condition wjj, = V sin a. The coefficient 
, of the principal diagonal are all equel and similarly in 
eaoh diagonal from upper left to lower right. Thus, for 
example, for X a 6, n"= 4; that is, = ± b/40, 

• 3"b/48, ± 5d/48, ± 7d/48 , the system "becomes:"" 



16.346i Tx - 14.345i T a - 4.250s l 3 - 2.284a r 4 = 2 °" ▼ sin a 

6.250s F,. + 16.346i r a - 14.346i T 3 - 4.250s ft = 2b V sin a 

4.284a Tj_ + 6.250s T a + 16.3461 T 3 - 14.346i r 4 = 2b 17 sin a 

3.4704 T x + 4.284a T a + 6.250s T 3 + 16.346i T 4 = 2b V Sin a 

with the solutions 

T x = 1.3344 t 7 sin a T 3 = '0.328 Q -t V sin a 

T s = C ,5f8 3 t V sin a 1^ = 0.179 x t V ain a 

The lift contributions from the four strips of the "in 4 ; 
are: Aj = 1.048, b t v cin a, A IX = 0.4785 b t V "in a, 

A IIX = 0.257s b t V sin a, A IT = 0.143, b t V sin a; 

the coefficients: c a = 3.76s sin cl. c w = 0.754 sin 8 a, 
Cjj = 0.92 3 sin a, and the distance of the center of pres- 
sure from tho leading edge s = 0.24, t. A comparison of 
the results "by this method and those "by- the Vortex-surface 
method will be made in the third section. 

Still simpler is the limiting case of the above sys- 
ten for n = 1. The conputation is then made with a vor*- 
texwith spanwiBe elliptic circulation distribution at the 
diatanco t/4 from the leading edge'. Tho single point 
consfdorod lios at x = 3/4 t. The center of pressure, ac- 
cording to the assumption thon, always-lies at s = t/4, / 
which is sufficiently accurate for rather large aspect ra- 
tios ?v > 3, as shown "by -the conputation with tho vortex: 
shoet . The circulation and the . coof f i ci ent s in this case 
can, in general, be explicitly expressod: 
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2h T sin a 



1 + £ J 1 + X s I 

TT 



tt X ain a „ _ 1 _ a _ _ 1 

c« = ■ ■ ' — — ■ ■ C_.i, — — r- Co , c r = — ca 

VTTT 2; 

With infinitely increasing aspect ratio there is thua oo- 

_a 

tainod a limiting value for c- : lin c n - ~- Bin a. From 

'X oo ** 2 

the potential thoory.ontho other hand, there is obtained 
for the win^ with infinite span c a = 2,-u sin a. This dif- 
ference is readily explainable "from the fact that elliptic 
circulation distribution (that is, factor 7, previously 
defined, equals n/4) has been assumed above. According to 
Betz, however, this factor for rectanfalar win^o "becomes 
larger with increasing \ up to the limiting value 1 for 
\ — * od . Hence, multiplying "by 4/tt, re have actually 

4 -lim / ^ 

tt X-*oD C a 3 ^ c a' po f 

It is now simple to refine this method by assuming a 
trigonometric series for r(y) s y = t? cos ©, r(<p) = T sin cp 
(1 + ai sin qp + . . . a n - sin n qp). Then, a^ain, equation (3) 
is obtained - thin tine for ore vortez with 6^ = l/X; 
that is, without the summation si^n. The dra? ic determined 

TT 

fron IT n p | /r(<p) w(qp) sin qp dqp, where w(qp) = 
tt % 

I dr(q?/ nn ^ there is obtained c_4 = 

2ttTjJ cos qp - cob qp t dqp w - 
o 

( 4 2 

+ - ai + - ai B + . . . . — j — ■ In tho examples there was 

^3 TT 7*"bt 

sot: r(tp) = T sin q> (1 + aj. sin qp) , cp x * = 30° ,cp a * = 60° , 

which corresponds to y* = b, y * s .2?. On figure 4 ■ 

4 a 4 

the lift-distribution coefficient is plotted aiainst the 
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aspect-ratio. It nay be seen that the experimentally ob- 
tained lift croef f icie'nts "(reference 8-) are very closely 

approached "by this simple aamputation. 



The following table is to be used in connection with 
figure 4 : 





°a 
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-f x 
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sin 8 
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. — J X sin 
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TT 

-4" 


4/7 


0.89 e 


sin 


a 


0.41 7 


sin 2 


A 


0.22 4 sin 


a> 


0.785 7 


4/3 


1-814 


sin 


a 


0 .78s 


sin a 


a 


0,454 sin 


a. 


0.788 e 


6 


4.202a 


sin 


a 


0.93a 

I 


sin 3 


a 


1.05i sin 


a 


0. 83G 7 



II. TORTUS SHE3T 



With the method of ciscrcte vortex filaments, a fur- 
ther refinement in tho lift distribution - hence an in- 
crease in the accuracy by incroasin*; the nunbor of vorti- 
ces — is practically excluded on account of tho computa- 
tion labor involved. It uas therefore larried to the Un- 
it n -» ib, ard an attempt was nade by analytical methods 
to restrict the computation work as far as possible On 
increasing n, the circulation contributed by each of the 
vortices and their distances apart become smaller up to 
the linitin-; case n — » &> . when the circulation distribu- 
tion bocor.es a surface distribution. The dimensions of 
fhls surface distribution Y(x,y) are circulation per 
unit chord, that is, centimeters per socond. An infinitesi- 
mal vortex "Y(x,y) dx inducos, according to equation (l), 
at any point of tho surfaco, tho • yolocity: 

+b/2 



-b/2 
+b/2 




I 
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whore a isaaaunod snail enough so that sin a « 0 and 

cos as 1. The equation would "be strictly true if the 
trailing vortices wero sbod in the diroction of tho win?. 
Into-; rating the first integral "by parts, thoro is a^ain 
obtained, on account of Y (i, y = =t "b/2 ) =0, ■ 

t- +Tj/2 

w A = - JL fax f j /(x*-x) a + (y^y*) _ + i I avCx^ a, 

&• - b / 2 

The condition wj^ = V sin a ^i*?es, after a transf ornation 
of coordinates, 

+i +i 

x =.t/2(i + e) r r [ /(£*-.t) a -fx a (n-Ti* ) a 

y = -b/2 ti "L' J L ce*-t) cn-*n*> 

- 1 -i 

+ -=-^- > ^li-JLl d t dTl = 4tt A V sin a (5) 
| oT] 

T^J s is r. two-dinensional linear integral equation of the 

first kind i a -rv — , On account of the sin i ralarity 

of the kernel at the t>oint considered for 5 £ , this 

Tl - 71* 

equation could not be solved eTen approximately since the 
usual nethod of approximating the kernel "by a polynomial, 
assumes continuity. Also the particular property of the 
kernel on which it depends l{- I*) and (T] -• T|*) could 

not be uaed for a solution nethod. 

The only recourse therefore is to sinplify the inte- 
gral equation "by au aerodynani cally reasonable assumption. 
It was thus assumed that the spanwise circulation distri- 
bution is the sane for nil £ . This results in a lowering 
of the order of the integral equation, although the flow 
condition can then no longer be satisfied ovor tho ontiro 
win*; but only on a straight line T)* = const. For obvious 
reasons 'elliptic distribution was aosunoi over the span, 
and T)* = 0; that le . tho flow condition was to do satis- 
fied on the center lino of the plate. For t^eee points of 
the surface then, the condition = V sin a is required. 
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and for other point's it Is expected'' that tho deviation 

fron this condition-la not- too lar^o* --' 

, . ■ 

Mtn Y(t.u) = v(e) ./i"T5 i i .aatii*a) = - 

. ol ' v 1 - u ■ 

anjL Tl* = 0, equation (5) simplifies to the following 
+ 1 +1 



d'Tl = 4tt X V sin a 

J 

-l -l 



In the Integration with respect to T|, on elliptic integral 

is a^ain obtained: 



+ i 



AC -if a " 



+ x " an = 



(t* -t) JT-f 



^» / A - "5 ~ sir." qpdcp, 



0 



with T] = cos cp, so thpt 



P r / U*-*) a +x a . C x n) + nlvrn d* - 

J 1 (f-t) v^TcP^TF 2 ■- 2 j 



2tt X 7 sin a (6) 



Corresponding Integral equations can also Ue set up for. 
airfoils with arbitrary plan form synmot ri cal with respect 
to the center line. If h( i) la the span, varying with 
the chord, equation (5) 'becoces: 

i 1 avct.-TV) i 4 , d71 _ 
tt^J— 5T~ Mil de dn " 



-i —i 

+ _ 1.- > ~ "H 1 . — v dt dTj = 4tt V sin a 

(Continued on p. 14) 
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Equation (5) is tbus reduced to a one-dlnenef onal integral 
equation. The kernel, it is true, has "become- more compli- 
cated, and the singularity for £ —* naturally, still 
remains. Bor the approximate solution of this equation, 
there is assumed for a series of functions which 

Birnoaum has used in his paper 

7(1)3 Ai JWi + Aa J*-? + + a < ffi 1 ? w 

The four undetermined coefficients A 1 ,A a , A„ A 4 are so 

determined that the integral equation is satisfied at four 
points £ x to £ 4 . Since tho kernel contains ^ , it is 

necessary to compute each tine for a deficit e aspect ratio. 
The integrations, on account of the cor.pli cated kernel , 
nust be carried out graphically or nunori cally . The pro- 
cedure of tie ver;- laborious computation thus, is the fol- 
lOTrin?: 

Che four basic functions , Y a , Y 3 , and Y 4 for a 

series of values of \ froo -1 to +1, are first comput- 
ed; then for thesame arguments for a definite A. the 
kornel functions k(X;|*,|) for the four points consid- 
ered. (In the examples carried ouf | x = - li = 0, 

£ 3 * = + i, U* = + 1; then, on account of 1: C- ) = 
k (£*,£) + it, k (-0.5, i) = k (+0.5, I) + tt.) Each of these 
kernel functions is multiplied by V x ([), Y a ((), Y 3 ([), 

■*«(£): 

(Continued from p. 13) 

For the elliptic win^ thero is obtained for elliptic circu- 
lation distribution over tho span with X ( £ ) ='b =a3C / 1 ^/l"^ i 

P ( z ass ziiig E ( * sm* e) 

- ; i l (T-i) V^ a (i-e a ) + (i*-i) a 2 

+ H ] _7liJ — d£ = 2rr V sin a; \ = 2sK 

2 JXvT-| a * 

This oqua^ion is solved approximately as in tho case of tho 
rectangular win^5. Tho limiting case A » *?ivos for 

1 +l 1 

2tt j f |*-£ Y( O&l = v ein a th. e solution of the potor.iic.1 

theory Y(£) = 2V sin a ~ with c a = 2tt sin a, c_ =» 
~ sin a, s = — t . 
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%1 


Kli - 


» 




i_ ^ t * t\ 
= k\l a it J 


Y f O 


ia 


(1) = 




i.Ct) 








CI) = 










*4 


(t) = 
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etc. to 

ii.(t) = kce/.t) 

Those 16 functiona must now "bo Integrated graphically or 
numerically : for exarple, "by theSInpaoa. rule: 

+i 

ik =j ' i k U) dt 

-1 

Since the k(|*,£) and hence the i(f) become infinite 
for £ — * 5 • it ^ i s ne coo cary in the quadrature to 8x- 
clucLe a region | -e<|<£*+€ and approximately de- 
termine "by anrlytical methods the Cnuciy principal valuo 



for J* i(|) dg. Tic thua finally :*avo tho coefficfonts 
I]^ for tho linear nonhono^onoous fly at or. of orations; 

4 

E^ Z 4^+U = 2tt X V sin cc. for \i = 0, 1, 2, 3 , . . . 

fron which the A V (X) can bo determined'. T7o then have: 

* / a 

Y ; 6 . Tl) - EuAy(X) V V (D Jl - Tl . in this nanner tho cir- 
culation diB^bri Dution was dotormincd for tho aspoct ratios 
A = 1/4, 1/2, and 1. 

Tor greater aspect ratios, A > 2, the 15 graphical 
quadratures are not required, but the approximate computa- 
tion can be carried out analytically. Por this purposo 
the elliptic integral must approximately evaluated for 

45° < q> < 90°. Since sin cp = X and U -£|< 2, 



16 
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cp = 45° Is tha smallest ar^iraent for X > 2. For this 
ran^e E(x, tt/2 ) was replaced "by E (x, tt/2) = 1 + 
0.44 Jl - x a . At first there was set E (z, tt/2) = 1 + 



^ -1^ ,/ 1 - x a ZV& V x v , but the computation then be- 
came so complicated that there was no adventa^e wined 
over the previous method; E is thus determined with an 
error which is approximately -3 percent for cp = 45°j +3.5 
percent (maximum value) for cp = 80°, and approaches zero 
as cp — » 90°. For cp = 90° , the first derivatives also 



dE fx 



a?rco: lim -"- - - » and 

x - i dx 

tt/s 

dE 



("■ f) 

lim =— lie / — ^dcp-»-oo 



lie f - z 
i - x iX x i J ,/i-p sin 3 cp 

For X > 2 the radical can "oo developed into a porcr S o- 
rio8 rhich converges for ell I , since - £| < 2. 

tfith 

e f X H^l ( * H N , = 

and 

1 + 0t44 _U!=il_ 



equation (6) then oeconep: 
+ l 



/ {^-r + s: (£ *" 5> "A «*- s;3 + nV ^>*t 



+i 



+ 



J (o.44 si?n (J* -0 + Y(£) dg = 2tt A v sin a 



- 1 

A*;ain substituting Y(£ ) from equation (7) and integrating, 

there ore obtained on the left side for the coefficients of 
A^, Aa i Aj, A« , the following functions of X and |". 



+ 0.280 j, j. (arc sin , 
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1 L16X B ' 32X 8 v 16 X B BX 3 ' 

+ ( 15 3\ ., a + / 15 _ 3 + _1\ |* 

^64X 8 16X 3 ' ■ \L28X B 16X 3 2 A' 

All 5 tt\ 

+ ( X + — _ ■ + — =■ + - ) 

>> 4A 64A 3 256A B 2' 
8 l.32X B ^64X B 16X 3 ' 

+ fx + i. - j- + —5—) r + a 

^ 4X 64X 3 256X 6 ' 4 
+0.140 06 (arc sin f + |* / 1 - I* 8 )] 

3 & L 12<3X* ^ 64 X s 128X 5 ' 

- (L + -L. - 1 - 0.09337 (l-r a ) 3/a l 

*2 16X 128X 3 2048X 5 ' * 

4 L128X* ^ 34X 128X" > 

- (h - -i. + -J: " \ r + n_ 

\2 16X 128A 3 2048 X 8 / 16 

+ 0.035 01 (arc sill £* - t* Jl - I* 8 (l - 2£* 8 )>] 

A^ain 8 definite value is takon for X and for four 
points (hero, too, the points chosen were | a * = 
la* - °' Is* = + if a n & l«* = + 1) thore are conputod 
tho 16 coefficient's of the oysten of oquati'ona 
Zu I,, (|li*) A, = 2tt X V sin a for W- 1. 2, 3, 4. fron 
which the coefficients- x A a* Ja.A, aro determined. 
In this nannor the A, for X = 6 and X = 2 aro conput- 



BACA Technical Memorandum ITo. 9 63 



od. For "K = 2, this computation, strictly speaking. Is 

at least for £ 4 * = 1, not valid "because -the oeri as into 

which fho radical was developed is no longer convergent in 
the Uniting case £ - 1. 

Finally, equation (6) waa'alao .cpnaldorod for the two 
linitinq cases A -» oo and A -* 0. ■ Bor tho win? of in- 
finite soan, -+ oo, E ( — ===== ? J "becomes oqual 

b +i 

n 

to 1 and tho equation goos-'ovor into j 7 ( t ) d£ = 

- 1 



2tt v sin a. Substitutinq the above expression for 

there is easily recoqnized as a solution *Y(£) = 2V sin a 

J ~~~ — ^ ' which is the distribution given "by the poten- 
tial theory. Since as A— *co, the spanwise distribution 
becomes Y ( T) ) = const, this solution sati sf ios -the flow 
condition at each point of the surface. The sane result 
la also obtained-when in the method of solution for \ > 2, 

the I\; • I*) are considered for very larsje X and this 

system is computed. Then there is also obtained Ax = 2, 
As = A 3 = A4 = 0. The lift coefficient rill then be c a = 

3 3 ^ 

-g- sin cc, and the moment coefficient Cm = sin a 

whereas, accordinq to the two-dimensional potential theory, 
c a = 2tt sin a and c m = sin a. This is a^ain due to 

the fact that elliptical epanriae distribution was assumed 
for tho rectangular win 1 ;; multiplying, subsequently, "by 
4/tt. the two results become identical. If an elliptical 
wins is considered and X is nade to approach infinity, 
thore is immediately obtained c a = 2tt sin a and c m = ^ 
sin a since the roforonce area for tho coefficients is 
it/4 b t. 

Of considerably groator difficulty is the limitin" 
case of the win? with infinito chord t — * ob, \ — * 0. 
Here the coordinntoa must "bo nade nondinensional through 

the span instead of through the chord as heretofore: x = 

^ v, x* =* ^ v*. Equation (6) with t -*■> oo then ?oes over 
into : 
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oo 

/{ 



Vi+(t*-t) 

Again it T7as sought to find for Y(v) a Horlos of func- 
tions with undetorninod coofficionts which would then be 
doterninod through satisfying the integral equation at sov- 
eral points v* . In this case, howovor, no sorios of func- 
tions could bo found which at tho leading edgo v — * 0, 
Increases as l//"r with -v —9 oo, corresponding to 
the solution for A = l/4 , which decroasoc approximately 
as l/v 3 . For this roason, only tio following single func- 
tions were investigated: r, (v) = —p= r and Y a (v) = 

J v + 3ir 

_ c 

• Since again tho integral could bo ovaluatod 

J v + Ev 5 

only graphically or nuneri cally (on account of 0<|v*-v|< oo 
a 3ories davolopcent of tho'-ornel could not be considorod), 
it was necessary, "before substituting, to anauno 3, and D, 
respectively , an firod and then, by quadraturoa, Set up 
A(3,v*) and C(E,v*), roopocti voly , for several valuos of 
v*. in order to linit tho integration intorval, it was 
a c ;ain necessary to np.ko another transf ornation : 

u = — i — . . u* = 1 

1 + v 1 + v* 




Since the integrand again "became 'singular at two points, t- 
nanely, nt u -*> u* and u — 1 (u = 1 corresponds to the 
leading edge, u = 0 to the trailing edge), the principal 
values hcd to be approximately detorcined by analytical 
methods. The ontiro laborious trial procous, howovor, not 
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with little success as tho effect of the various coeffi- 
cients on the result was too difficult to estimate. As an 
approximation, it is possible to sot at nost Y(v) ~ 
1 12 

~ - ; for the ran*?e 0.1 < v < od, the downwash 

•fr + 2.4 v 3 

orror with respect to v sin a ■chen amounts to about 7 
percent - this error, however, strongly increasing toward 
tho loading ed^o (v < 0.1 ). Tho lift would then amount 
to A = 0.789 b a V a sin a, r.nd the contor of uroosuro 

rould lie fit s = 0.219 b. 



III. RESULTS 



3-" the methods deacribod, the chordwise lift distri- 
bution was conputcd for elliptic soanwise distribution for 
five aspect ratios, namely, A = 1^4 , l/2 f 1, 2, and 6, and 
graphically interpolated for arbitrary X. The results 
are presented in figures 5 to 10 and are tabulated in the 
appendix. In figure 5 the coefficients for the circula- 
tion functions are plotted against A. A x Increases mon- 

otonically and for lar«;e aspect ratio approaches 2 as the 

asymptotic value. The absolute values of the other coef- 
ficients increase up to a maximum at about A = ■g, then 
drop rapidly to zero; A, and A, are always positivo, 
and Aq and A 4 , negative. Tho smaller A is, the less 
rapidly do tbe A^ converge, sothat to obtain the same 
accuracy as for larfife aspect ratios, a longer sorias of 
functions for the circulation must be assumed. Tho curve 

A x /bA shows the increase in the circulation in tho neigh- 
borhood of the leading od%e for constant span as 3 func- 
tion of X, since lir Y(x)=lLi I r ~7= • here, too, 

x - o v X y x 

the maximum lies between A = J and X = 1. The valuo 

for A = 0 obviously is in orror, fror. which fact it may 
ho soen that the ^iven approximation for Y(y) does not 
correctly represent the behavior in tho neighborhood of 
the loading ed^o. If tho chord is held fixed a:id tho span 
varied, Ai itsolf -sives the increaso since lim = 

, x -*■ o 

A 1( /t/x. 



On figures 6 and 7 the circulation distribution 
*Y(£)/V sin a and tho pressure difference hotweon the 
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lower and upper sidos of tho plate referred to the dynamic 



chord. In figure 6 the abscissa refers to the chord, and 
ir. figure 7, to the span. Ir. tie first representation the 
limiting case "K = 0 coincide3 with the coordinate axes, 
so that tho lift distributions for on" aspect ratio lie3 
between tho axe3 and the Uniting case X = co. In tie 
second representation tte Uniting case A = o» coincides 
with axis of ordinates. This method of plotting is par- 
ticularly susceptible to error3 in thi circulation distri- 
bution and show3 thnt tt A = J, a snail error is to be 
assumed through inaccuracy of one of the graphical quad- 
ratures or the approximated principal value. Similarl;^ , 
the limiting case X = 0 appears as only a ver^ - rou'Sh ap- 
proximation since intersection of the curves with each 
other is very improbable. It is seen, howovcr. that for 
very small aspect ratior a further lacrosse in tio chord 
has only a email effect on the circulation distribution, 
either in tho neighborhood of tho loading ei^e or - owir^ 
to tie strong decrease — farther toward ths rear. This is 
aeon osrecially clearly from tho curve for A/p b a V 3 sin 
(fi?. 8 v , which shows how the tot".l lift increases when 
■cho spin is held constant in? tho chord is varied. 

On figures 0 and 1C the lift, dr? , nnd moment cooff- 
cients are plottod as function3 of A. Thoy are computed 
from the valuos of Aj, A a , — 3 , A 4 a3 follows: The lift 
according to tho Kutta-Jourowslcy theorim is: 



pres sure 



Pu ~ ?o 



— ar'e plotted against the 



— 7 sin a 



+b/2 t 



A = 



P 




-b/2 



0 





Similarly, the 



moment about 



tho leading od-~o is 




+b/2 

r 



t 

n 



-b/2 . o , 




finally the in- 
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duced dra"! is 

+"b/ 2 * 



with 



*i = P J* dy J Y(x,y) w(y) dx 



-V 2 o 

where 



dw(r) = i / dx — ? 

4tt J ay y _y 

-b/2 



(x.y) -| ' ' ^ 8 



yi - ^ A V Y v (x) 

re nave 

4 

inisi = - 1 E„ a, Y,.(x)«= const 

- a x 2d w 



on account of tie olliptic epr\nwise distribution, so that 
v x = 21 1^ + | A a + | A 4 j -Mid 7 4 = 2| P t a (a x + I A a + 

i A 4 j M ft ionco e wl = E_ -i^ (^A x + | A a + | A 4 ) . 

c _a 

Eonce in this case also, = tt X . Cn account of t]?e 

°wi 

elliptic circulation distribution, the factor v = tt/4; 
hence multiplying by — ) to tako account of the v 

TT 

factor riich incroasos with A, (dc a /dx) Q incroasea up to 
2tt as \ -» oo. (The factor v was takon fron the dis- 
sertation "by A. Pets.) The induced dra- - ; coefficient in- 
croaEes with increasing A, has a r.axir.un at about "K = 2 
and then, since the dra? rer.ains finite vhile t>e area "be- 
couos larger, drops to zero. The position of the center 
of pressure is obtained fron s = (0,3/ c a ) t. This curve 

rapidly approaches tho asymptote. At X = 3, the devia- 
tior fron the Uniting value s = 0.25 t is only 5 per- 
cent (fi^. 10). 

The a<jreonent of tie conputation results "by the vortex- 
filament r.ethod aith those "by the vortex-sheet rethod is 
surprisinqlp ^ood. Even with tso vortex filar.ente tho de- 
viations in the coefficients ar» snail, rhoreas with four 
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vortex filaments the deviations become lar«e .only for very 
deop. .plates with X< 1. (For X = 1: Ac a /c n « - 0.4 

percent, Ao n /c in »+ 3.4 perrerit; for' \ = 1/-2 : Ac a /e a = 

- 1,6 percent, ACju/cn s + 90 percent.) The individual 

lift portions contributed "by the four strips of tho sur- 
face do not agree so well; for example : 

Four vortex f ijLaments Vortex eheot 

Lift fron 0 - 3t/l6: 0.711 "b t Ysina 0.665 h t V sin a 

3t/l6 - 7t/l6 : 0, 2 64 n 0.274 11 

7t/l6 - llt/lo: 0.141 » c.144 11 

llt/16 - t : 0.071 " 0.378 11 

i 

Tho circulation of the foremost vortex alrays cones out too 
hi^h, and that of the other vorticos too low. In obtain- 
ing the moment this error is partially compensated oy tho 
consi deration that too 2.ar°;e lever arnc r.re used for the 
three rear vortex f llanont s , xrhi ch do not lie at the cen- 
ters of gravity of Ajj, Ajjj, Ajy 

» 

With both nethodo the assumption of elliptic distri- 
butiea over the span - vhich assumption makes possible the 
solution of the integral equation in tho case of tho vortex 
sheet - should be the greatest source of error. For this 
reason, too low lift coefficients ore also obtained. Tho 
subsequent multiplication "by the factor v does not ap- 
pear to help sufficiently. ."-icordin? to tho pressure dis- 
tribution . measurements of H. winter rhich, howovor , are od- 
tainod for the square plate only, the distribution over the 
span at the loading ed^e is approximately olliptic, but 
farther toward tho rear - alrostup to the ed-^e - It is 
constant, tho ed<;e di sturbanees which ari se from the sharp 
ed^es of the investigated plate, horever, not bein*5 taken 
into acoount. 

, APPENDIX 



I. 'Several vortex filaments: Tto vortex filament8 
with elliptic circulation distribution over the span at 
x = t/8 and x = 5t/8. 
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X = 4 Ai =0.3571 T)t7 sin a 

A x 2=0.0277,3 l)tV sin a 

X = 1 Aj =0.6132 T)t7 sin a 

AjjsO.1050 bt7 sir. a 

A = 2 Aj =0.9410 bt7 sin a 

At I= 0.2434 "bt7 sin a 

A = 6 Aj =1.4280 l>t7 sin a 

AjjsO.4557 T>t7 sin a 



o a =0.770 sin a 



0,3=0,124 



sic a 



c a =1.435 sir a 

CjjjsO.285 sin a 

0^=2.269 sin a 

03=0.540 sin a 



c a =3.767 sin a 



CjasO.923 sin a 



8=0.161 t 



a=0.198 t 



8=0.228 t 



s=0.246 t 



Tour vortex filaDenta with elliptic circulation dis- 
tribution over the g-pen r.t x = t/lo , 5t/l6, 9t/l5, r>nd 
13t/l6. 















Vortex shoot 








= * 


±1 


=0.5080 bt7 


sir 


a 


(0 .2979) 


o a =0.772 


sin a 






A II 


=0.0557 " 






(0.0724) 


0^=0.379 


siri 3 a 






A III 


=0.0165 " 






(0.0154) 


c n =0.101 


sin a 






Aiy 


=0.0058 » 






(0.0067) 


s =0.131 


t 


A 


= 1 


A I 


=0.4876 "bt7 


sin 


a 


(0.4838) 


c a =1.441 


sin a 






A II 


=0.1433 n 






(0.1470) 


c w =0.551 


sin 3 a 






A III 


=0.0623 11 






(0.0627) 


c n =0.265 


sin a 






A IT 


=0.0273 ■ 






(0.0299) 


s 40.184 


t 




= 2 


A I 


=0.7111 t)t7 


sin 


a 


(0 . 6852) 


c a =2.374 


sir. a 






AI I 


=0.2544 » 






(0.2739) 


0^=0.897 


sin 8 a 








=0.1405 11 






(0.1442) 


032=0.528 


sin a 






A IV 


=0.0710 11 






(0.0782) 


s =0.222 


t 
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Vortex sheet 

A = 6 Aj =1.0480 btT sin a ('1.0068) e a =3.770 sin a 

An =0-4385 H (0.4425) . <^r= Q ' 754 Bln a a 

A III= 0 »2576 « (0.2525) e n =0.,923 sin a 

Ajy =0.1406 M (0.14 65) a =0.245 t- 

J I Tort.ftT -.- surface : ExanploB for a < 2 
a = i J x 6.3138 J a =0.6725 Jfe = -0.7284 J 4 =0.5566 

J B = 8.2 5 2 3 J 6 =2.4674 J 7 = -1.1575 J B =0.6169 
J 0 =9.4427 J 10 sJ. 0623 J a x = -0.7284 J la =0.6772 
J 13 =10.0592 J 14 =5.2057 J lB = 0.2 5 8 1 J 16 =1.4505 
Ai = 0.3435 v sin ct 0^=0.3839 sir, a A x /J\=Q .687 V sin a 
A, =-0.3307 " 0^=0. 1"25 cin 8 a total lift: 

A, = 0.2994 " c n =0.0287 sin a A = 0.7778 p b a V 

A, = -0.1542 11 s =0.0739 t 

J x = 6.8189 J a =0.6 345 J 3 = -1.0013 J 4 =0.6956 
J B = 8.6 5 0 8 J B =2.4574 J 7 = -1.5644 J 8 =0.615? 
J e =9.7617 J 10 =4.2703 J 12 = -1-0013 J ia =0.5381 
J 13 =10.38i6" J 14 =5.8817 J 1B = 0.5725 J ie -1.7546 

i-i = 0.5793 v sin a c a = 0.7846 ein a A x ^^=0.819 V Din rt 

A a = -0.4802 «' c w =0.3919 ain^a A = 0.7846 p b 8 V 8 

A 3 = 0.3710 n 0^=0.0922 sin a 

A, = - 0 . 1 6 9 7 n s -Q.H75 t 
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X = 1 J x = 8.1785 J a =0.1265 Jj * -1.1401 J 4 =0.8755 

S B = 9.6696 J 6 =2.4674 J 7 = -1.98 64 J 8 -0.6169 

J 9 =10.6560 J 10 =4.8083 -1.1401 J ia =C.3583 

J 13 =11.3206 J 14 =7.0036 Ji B = 1.3577 J ia =2.4958 

A x = 0.8182 V sin a c n =1.4456 sin a Ai^T=0.818 V sin cc 

A a = -0.4424 " c v =0.6562 sin 8 a A = 0.7233 p b 3 7 2 

A 3 = 0.2440 ■ o n =C.2563 sin a 

A 4 =-0.0852 » s =0.1771 t 

Exanples for X > 2 

*«2 F x =3.6683" F 3 = -0.4078TT F 3 = -0.5899tt F 4 = 0.4204TT „ 

F B =3.9707it F 8 = 0.7854tt F 7 = -1.1236tt Fb = 0.1954tt 

F p =4.1905tt F 1C = 1.9786tt 7 lt = -0.5898tt F 13 = -0.0277- 
F i3 =4.3408" F i4 = 3.1208tt F 15 = 0.9732ttF 16 = 1.2815tt 

Ai = 1.0790 V sin <r. c a =2. , 3629 sin a A 1 / v /^=0.763 V sin a 

A a = -0.2389 " 0^=0.8886 sin 8 a A = 0.5937 p l) 8 T° 

As= 0.0834 " c m =0.5288 sin a 

A 4 = -0.0149 » s =0.2238 t 

*=6 J 1 ! =7.6668tt F a = -2.3693" F 3 = -1.5710" F 4 = 0.9304tt 

F 6 =7.8924" F 6 = 0.7854" F 7 = -3.1037" F a = 0.1964" 

F 9 =8.0425" F 10 = 3.9401" Fj 1 = -1.5710H J ig = ~0.5377tt 

■ F 13 =8.1334" F 14 = 7.0466" F 15 = 2.9898" Fi 6 = 3.2517tt 
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A x = 1.5449 v sin % c a =3;6946 Bin a A. 1 /fT^Q.6Zl V sin a 

A,= -0.09344 » e w =0.7242 ain s a A =0.3079 p if V 3 

A 3 = 0.0380b " c m =0.9002 sin a 

A 4 = -0.0063 o 11 s =0.2437- t 

X = co A 1 =Z7 sin a Cq> = -g- gin a instead of c a = 2tt sin a 

A a =A 3 =A 4 =Q o m = Hi Bin a 11 o n = j ain a 

s = 0.25 t 



Tr-nslalion "by S. Reiss, 
ITr.ti 0 n. n .I Advisory Committee 
for Aeronaut! ca . 
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Jigure 4.- Lift distribution ai a function 

of A for a vortex filament, 
a, r (^rsin tfj rOPj-rainVU+aj, sin^). 
c,te«t result8 of Winter. 
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Figure 7,- Circulation dietribution 
over the chord; b= conet . 



Figure 3.- Circulation dietribution 

of the forward three 
vortex filaments at x=t/l6, 5t/l6, 
9t/16;X= l. 
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Figure 5.- Coefficients cff 

circulation 
functions for various aspect 
ratios . 
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Figure P.- Circulation distri- 
bution over the 
chord; t = const, 
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Figure 8.- Total lift as a 

function of the chori. 
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A = l/t 

Figure 10,- Moment coefficient and distance of center of 

pressure from leading edge as functions of X, 
(points x computed with four vortex lines). 
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